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Abstract. In this paper, we investigate power allocation that maximizes 
the secrecy rate of orthogonal frequency division multiplexing (OFDM) 
systems under arbitrarily distributed inputs. Considering commonly as- 
sumed Gaussian inputs are unrealistic, we focus on secrecy systems with 
more practical discrete distributed inputs, such as PSK, QAM, etc. While 
the secrecy rate achieved by Gaussian distributed inputs is concave with 
respect to the transmit power, we have found and rigorously proved that 
the secrecy rate is non-concave under any discrete inputs. Hence, tra- 
ditional convex optimization methods are not applicable any more. To 
address this non-concave power allocation problem, we propose an effi- 
cient algorithm. Its gap from optimality vanishes asymptotically at the 
rate of 0(1/a/A), and its complexity grows in the order of 0{N), where 
N is the number of sub-carriers. Numerical results are provided to illus- 
trate the efficacy of the proposed algorithm. 
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1 Introduction 

In recent years, many privacy sensitive wireless services, such as pushmail, mo- 
bile wallet, Microblogging, etc, have become more and more popular. But, due to 
the broadcast nature of wireless channels, security problems and challenges are 
also accompanying with the growing up of those privacy services. The security 
of wireless communications is commonly supported by cryptographic techniques 
employed at upper layer. However, this traditional method faces several chal- 
lenges, such as the emergence of new cracking algorithms and increasing compu- 
tational capability of eavesdroppers. Recently, physical layer security, a method 
that can supplement upper layer security, has received considerable attentions 
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Physical layer security was firstly studied from an information-theoretic per- 
spective in [2 , where the concept of "wire-tap channel" was introduced to il- 
lustrate the channel with three terminals, transmitter, legitimate receiver and 
eavesdroppers, and secrecy rate was defined as the maximum achievable data 
rate from the transmitter to its legitimate receiver while keeping the eavesdrop- 
per completely ignorant of the secret massage. Later, the research in this field 
was extended to various scenes, such as Gaussian wire-tap channel [3]- [4], multi- 
ple input multiple output (MIMO) channel [5]-[2, orthogonal frequency division 
multiplexing (OFDM) channel [8 -[iP, etc. 

Recently, OFDM-based secure communications have obtained much atten- 
tion for its capability of countermining the dispersive of wideband wireless chan- 
nels and enhance secrecy rate f8l-[TT]. Optimal power allocation of OFDM-based 
wire-tap channels is investigated in [8]- [10] under Gaussian inputs. In practical 
system, due to Gaussian inputs' infinite peak- to- average ratio, finite discrete 
constellation^ such as PSK, QAM (see Fig(T](a)), are used instead. In this pa- 
per, we investigate optimal power allocation for OFDM-based wire-tap channels 
under arbitrarily distributed channel inputs. 

We have found and rigorously proved that the secrecy rate is non- concave 
under any discrete constellations, while the secrecy rate achieved by Gaussian 
distributed channel inputs was found to be concave with respect to the transmit 
power [8]- [To]. Therefore, the optimal power allocation strategy for OFDM-based 
wire-tap channels with Gaussian inputs 0-[IO] is not applicable any more to the 
considered problem. To address this non-concave power allocation problem, we 
propose an efficient power allocation algorithm. Its gap from optimality vanishes 
asymptotically at the rate of 0(1/V^), and its complexity grows in the order 
of 0{N)^ where N is the number of sub-carriers. Numerical results are provided 
to illustrate the efficiency of the proposed algorithm. 

The remainder of this paper is organized as follows. Section [2] presents the 
system model and power allocation problem is formulated. Optimal power allo- 
cation for arbitrarily distributed inputs is given in Section [S] Numerical results 
and conclusions are provided in Section S) and Section [S] respectively. 

2 System Model and Problem Formulation 

Consider an OFDM-based wire-tap channel with a transmitter, a legitimate re- 
ceiver and an eavesdropper, where the eavesdropper intends to extract the con- 
fidential message transmitted from transmitter to legitimate receiver (see Fig. [1] 
(b)). There are N sub-carriers and the transmitter's signal in each sub-carrier 
follows an arbitrary but predetermined distribution, which can be either continu- 
ous constellations, such as Gaussian distribution, or finite discrete constellations, 
including PSK, QAM, etc (see Fig.[I](a)). 

The transmitted signal over the ith sub-carrier is denoted as Xi , described as 

^ The words "distribution" and "constellation" are used alternatively throughout the 
paper. 
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Fig. 1. (a). Finite discrete constellations: QPSK and 16QAM, (b). OFDM-based wire- 
tap channel 



Xi = ^/PiSi.i = 1, . . . , A^, (1) 

where pi is the power ratio between transmission signal Xi and the noise, and Si 
represents the normalized channel inputs with predetermined distribution. Then 
power constraint can be readily shown to be 

i=l 

where P is total available transmit power. 

The received signals at the legitimate receiver and eavesdropper are given by 

yi = hi^iSi^Wi, z = !,•••, A/', (3) 
Zi = Qi^/PiSi + Vi, i = 1, • • • , A^, (4) 

respectively, where the Wi and Vi are zero-mean complex Gaussian noises with 
unit variance; hi and gi are the complex channel coefficients of ith sub-carrier. 
According to the information theoretical studies of [8^, the secrecy rate from 
transmitter to its legitimate receiver is 

N 

^[I{si] hiy/plsi + Wi) - I{si]gi^iSi + + , (5) 

2=1 
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where [x]~^ = max{x,0}, and I{x;y) denotes the mutual information between 
random variables x and y. The expression in (j5j) is quite illuminating: the secrecy 
rate of each sub-channel is non-negetive; if it is positive, it is exactly the data 
rate difference of the legitimate and eavesdropping channels; the total secrecy 
rate is simply the sum secrecy rate of all the N sub-carriers. 

For fixed constellations of {si}fLi, we need to optimize the power allocation 
to obtain the maximal secrecy rate. The optimization problem can be cast as 
follows, 

N 

R* = max Rs{p) = ^ Y.[I{si; hi^Si ^ Wi) - I{si;gi^Si + + 
p > 

where p G TZ^ is transmit power vector for N sub-carriers, i.e., p = {pi,P2, •••7 Pat}, 
and denotes the optimal value. For the facility of the following analysis, let 
Rs,i{Pi) denote the ingredients of Rs{p)^ i.e., 

RsAPi) - [I{si]hi^iSi^Wi) - I{si;gi^iSi^Vi)]^. (7) 



3 Optimal Power Allocation under Arbitrarily Distributed 
Channel Inputs 

3.1 Non-concavity of the Secrecy Rate Rs{p) 

If Si follows Gaussian distribution, the secrecy rate Rs{p) in (j6j) has explicit 
expression [4 , i.e., 

1 ^ 

It can be checked that Rf{p) is a concave function of p. Hence problem (j6j) is 
a convex optimization problem. The ingredients of Rf{p) are illustrated in the 
left part of Fig. [2l One can observe that log2(l + l^ipPi), log2(l + l^ipPi) and 
log2(l + \hi\'^Pi) — log2(l + Idil'^Pi) are all concave, provided that \hi\'^ > 
In [8]-[T0^, utilizing the convexity structure of problem ([6]), authors obtained the 
optimal power allocation shown as follows, 

{ -(|/^,P + |^,P) + v/(|/l,P + |^,P)2-4|/^,|2|^,|2 -+|g^l^-| ^ 
2\hi\^\gi\^ 
,if \h,\' - \g,\' > U (9) 
, others, 

where the Lagrange multiplier u is chosen to meet the power constraint. 
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One may expect that the concavity of Rs{p) still holds under general input 
distributions. Unfortunately, our investigation shows that this is not true, which 
is formally presented in the following proposition. 

Proposition 1. The secrecy rate function Rs{p) under any finite discrete con- 
stellations is non-concave with respect to p. 

Proof. When pi = 0, one can derive I{si;yi) = I{si]Zi) = 0; when pi = +oo, 
we have I{si;yi) = I{si;Zi) = H{si)^ where H{x) is entropy of x. Therefore, 

According to [15] , 

dl(s: ./ps + n) , , , , 

^ = MMSE(p). (11) 

MMSE(p) is defined as 



MMSE(p) = E[\s - E{s\^s + n)p]. 



(12) 




where E[x] is the expectation of random variable x; ]E[x|?/] is the conditional 
expectation of x for given y. By (fTTj) , the derivative of Rs i{pi) at = is given 

bil 

i?;,,(Pi)U=o= [|/ii|'-|5in^>0, (13) 

which indicates that there must exist a > that Rs,i{Pi) > 0. According to 
the Lagrange's mean value theorem [17 , it must have a point pi G [pi, +oo] with 
negative slop Rg^^{pi) < 0. 

Assume Rs^i{pi) is concave, then the inequality Rs,i{pi) < Rs,i{pi)-\-Rs i{Pi){Pi~ 
Pi) holds [18j, which indicates Rs^i{-\-oo) = — oo. This contradicts Rs^i{-\-oo) = 0. 
Therefore, the concavity assumption is not true, and Proposition 1 holds. 

Two evidentiary examples are provided to illustrate Proposition 1. 

The first example is QPSK. The curves of /(s^;^^), I{si]Zi) and Rs,i{Pi) 
versus pi are shown in right part of Fig. [21 and they are in accordance with the 
statements in the proof of Proposition 1. 

The second example considers a 4 points PAM constellation with non-uniform 
spacing. The probability mass function of Si is given by 

^ Only the sub-carriers that satisfy \hi\'^ > |^^|^ are considered, as Rs,i{pi) = for 
those sub-carriers with < which do not affect the concavity of Rs{p). 
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where L is a normalization parameter to maintain unit variance. Figure [3] shows 
the secrecy rate Rs,i{Pi) for this case. It is interesting that the Rs,i{Pi) has two 
peaks. Hence, it is definitely non-concave. We note that the mutual information 
I{si] Hi) and /(s^; Zi) are concave with respect to Pi in linear scale [16]. 

3.2 Optimal Power Allocation Solution of Problem ([6]) 

Although problem (6) is non-convex, there are still some efficient algorithms to 
solve it and obtain near-optimal solutions. One of them is the Lagrangian duality 
method [18]. Some recent studies [12] -[14] showed that, for certain non-convex 
structures, asymptotic optimal performance can be achieved by this method. 
The Lagrangian of problem (|6]) is given by 

1 ^ ( I ^ \ 

L{p,u) = — ^[I{si;hi^/p~iSi -^Wi) - I{si;gi^iSi^Vi)]^ ^^[^ ~ N ) ' 

"'(15) 

where u is Lagrangian dual variable. The corresponding dual function can then 
be written as 

g{u) = max L{p,u). (i6) 
Hence the dual problem formulation of problem (|6]) can be readily expressed as 

= min g{u), (17) 

n>0 ^ ^ 

where denotes the optimal dual value. Since the objective function of primal 
problem (|6]) is non-concave, there is a positive gap between i?* and I)*, i.e., 
D* — R* > 0. However, according to the recent studies of Luo and Zhang [12], 
[13], asymptotic strong duality holds for problem (j6j), i.e. the duality gap D*—R* 
goes to zero as ^ 00, as is expressed in the following proposition: 

Proposition 2. // the channel coefficients gi and hi are Lipschitz continuous 
and bounded in the sense 

\M - \hj\\ < ij-^^y i,j e {1,2, ...,N} (18) 

\\9i\-\gj\\<Lj-^,\fi,je{l,2,...,N} (19) 
where Lh,Lg > is the Lipschitz constant. Then we have 

0<D* -R* <0 f j • (20) 
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Proof. According to (pT]) and ([12]), we have [15] 



< + = MMSE(p) < E[|s|2] = 1, (21) 



which imphes that 

^^^''^^^ - \h,\'mSE{\hi\^p) - |5,pMMSE(|5,|2p) < \hi\\ (22) 



dp 

= 2p\h,\mSE{\h,\'p) < 2p\h,l (23) 
= 2p\g,\mSE{\g,\'p) < 2p\g,l (24) 

d\gi\ 

where Rs^i{p) is defined in the expression of ([7]). 

According to equation (|22]) and the Lagrange's mean value theorem [17^, we 
have 

|i?,,,(p)-i?,,,(pOl < M^\\P-P Woo- (25) 

On the other hand, combining ([23]), ([2l]), ([21]), ([H]) and ([19]), using chain rule 
and the Lagrange's mean value theorem ^17^, we have 

\R,4p) - Rs,j{p)\ < {2p\hi\Lh + 2p\gi\Lg) (26) 

As \hi\^ \gi\ and p are bounded, according to (|25]) and (|26|) . there must exist an 
L such that 



|i?,,,(p)-i?,,,(p')| <L(^^ + |b-p'|U) (27) 



i.e., the secrecy rate Rs,i{p) is Lipschitz continuous. Hence according to Theorem 
2 of ^2 5 duality gap between and i?* vanishes asymptotically in the order 
of 1 / >/]V, which is expressed as 

< - R* < O(^), (28) 
~ WN 



and Proposition 2 holds. 

The procedures to solve ([16]) and ([T7)) are provided in the following. 
For each fixed problem ([T6]) can be decoupled into N independent sub- 
carrier problems 

g{u) = max L{pi,u), 

Pi>0 
N 

= ^max[[I{si;yi) - I{si]Zi)]^ - upi] ^uP. (29) 

2=1 
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Table 1. 

Algorithm : Lagrangian dual optimization method 

Initialize u 

repeat 

for i=l to N 

find Pi = argmax[[/(si; ^i) - /(s^; Zi)]"^ - upi] + uP. 
Pi 

end 

update u using bisection method. 
until u converges. 



While the sub-carrier problem in (|29|) is still non-convex, it has only one variable 
Pi and can be solved by simple one dimension line search. As the dual function 
g{u) is convex in its subgradient g\u) = P — "iZiLiPi^ where p* is optimal 
solution for problem (p!6|) with fixed is an increasing function in u. Hence 
bisection method can be used to solve dual problem ([T7j), so that either u = 0, 

p* is satisfied. TableiU summarizes the 

algorithm. 

The complexity of this algorithm is A^^log2(^), where is the accuracy 
of one dimension exhausitive search to solve (fT6|) and Cd is the accuracy of the 
bisection search to solve (pT|) . Since its complexity is linear with respect to the 
number of sub-carriers N ^ it is quite convenient for practical large values of N ^ 
such as 64^4096. We note that the complexity of solve ([6]) exhaustively is 
which is exponential in N and thus unrealistic. 

4 Numerical Results 

In this section, we provide some simulation results to illustrate the performance 
of our proposed power allocation algorithm and show how different channel input 
distributions affect the secrecy rate and power allocation results. 

We first consider an OFDM-based secure system with N = 128 sub-carriers. 
The secrecy rate versus total power constraint under different power allocation 
strategies and channel input distributions are illustrated in Fig. IH Two reference 
strategies are considered to compare with our strategy: the optimal strategy 
under Gaussian inputs, i.e., ([9]), which is denoted by "PA of (|9])" in Fig. Hj the 
equal power allocation strategy, which equally allocates total power among the 
sub-carriers that satisfy \hi\'^ > \gi\'^ and is denoted by "equal PA". 

Seen from Fig. IH higher secrecy rate can be achieved for QPSK and 16QAM 
by our proposed optimal power allocation strategy, especially when the power 
constraint P is quite large. However, equal power allocation and power allocation 
of ([9]) can be quite bad under finite discrete constellations, and the secrecy rate 
even drops to zero for large value of P. Actually, when P is large, secrecy rate 
under Gaussian inputs can be approximated by 
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Fig. 4. The secrecy rate versus total power P 



30 



(P) = ^ log2(l + \h^?P^) - log2(l + \9^?P^) ] + 

(30) 

which is independent with power allocation pi. So equal power allocation works 
asymptotically optimal under Gaussian inputs when P is large, as shown in Fig. 

n 

The power allocation solution of the proposed algorithm is shown in Fig. [5] 
and Fig. [6l respectively, under QPSK and Gaussian inputs with = 4. When 
the power constraint P is small, most transmit power is allocated to the stronger 
sub-channels, the channels with larger \hi\^ — \gi\^ (Channel 2 and Channel 4 
in our simulation example). However, as P grows, the transmit power allocated 
to the weak sub-channels grows. Under QPSK input signals, the transmit power 
allocated to every sub-channel should stop increasing when P is very large. But 
the transmit power under Guassian input signals keeps increasing. 



5 Conclusion 



In this paper, we have obtained the optimal power allocation for OFDM-based 
wire-tap channels with arbitrarily distributed inputs. While the secrecy rate 
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Fig. 5. Power allocation results versus P under QPSK inputs. 
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Fig. 6. Power allocation results versus P under Gaussian inputs. 



achieved by Gaussian distributed channel inputs is concave with respect to the 
transmit power, we have found and rigorously proved that the secrecy rate is 
non-concave under any practical finite discrete constellations. A power allocation 
algorithm has been proposed, which is asymptotic optimal as the number of 
sub-carrier increases. Our numerical results show that more transmit power may 
results in a huge loss in secrecy rate, which is rarely seen in previous power 
allocation studies. This indicates that optimal power allocation is quite essential 
in practical studies of physical layer security. 
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